Superconducting Vortex with Antiferromagnetic Core 
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We show that a superconducting vortex in underdoped high Tc superconductors could have an 
antiferromagnetic core. This type of vortex configuration arises as a topological solution in the 
recently constructed SO(5) nonlinear a model and in Landau Ginzburg theory with competing 
antiferromagnetic and superconducting order parameters. Experimental detection of this type of 
vortex by /iSR and neutron scattering is proposed. 
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One of the most striking properties of high Tc super- 
conductivity is the close proximity between the antifer- 
romagnetic (AF) and the superconducting (SC) phases. 
While there are a number of theories linking the mi- 
croscopic origin of high Tc SC to antiferromagnetic corre- 
lations, it is natural to ask if the close proximity between 
these two phases could have any macroscopic manifesta- 
tions. Recently, a unified theory of AF and d-wave SC 
in the cuprates has been constructed based on an SO (5) 
symmetry. In this theory, the AF and the d-wave SC or- 
der parameters are unified into a five dimensional vector 
(ni, n2, ns, n4, ns) called a superspin. The AF order pa- 
rameters Na correspond to the {n2,n3,n4) components, 
while the real and imaginary parts of the SC order pa- 
rameter A correspond to the (ni, ns) components. It was 
shown that the chemical potential induces a first order 
superspin-fiop transition where the superspin abruptly 
changes direction from AF to SC. 

The SO (5) theory predicts a spin triplet pseudo Gold- 
stone boson associated with the spontaneous breaking of 
S0(5) symmetry in the SC phase these can be iden- 
tified with the recently observed resonant neutron scat- 
tering peaks in superconducting YBCO 0. Physically, 
these modes corresponds to Gaussian fluctuations of the 
superspin. However, it was noted |^ that the S0(5) 
theory also admits a special class f topological solutions 
called meron configurations. In this configuration, the 
superspin lies inside the SC plane far away from the ori- 
gin, and the SC phase winds around the origin by 27r. As 
the origin is approached from the radial direction, the 
superspin lifts up from the SC plane into the AF sphere 
in order to minimize the energy cost of winding the SC 
phase. The result is a SC vortex with an AF core 

The existence of SC vortices with AF cores leads to 
non-trivial macroscopic consequence which we shall ex- 
plore in this paper. We first present detailed analytic 
and numerical solutions of the SO (5) vortex. We also 
study a more general Landau- Ginzburg (LG) theory ob- 
tained from the SO (5) theory by relaxing the constraint 
on the magnitude of the superspin. This theory describes 



AF and SC order parameters in competition with each 
other. We show that even if the SC state wins in the bulk, 
under appropriate conditions a nonvanishing AF compo- 
nent can occur inside a SC vortex core. The nature of 
the condition leads us to conclude that a SC vortex with 
AF core should only be realized in underdoped high Tc 
superconductors, not in the overdoped ones. We believe 
that the nature of the vortex core has nontrivial implica- 
tion for the physics of high Tc superconductors in a high 
magnetic field. In recent experiments, Boebinger et al. 
Q find that insulating and normal phases appear upon 
destruction of SC by a high magnetic field in underdoped 
and overdoped materials, respectively. This observation 
could be intimately related to the insulating/normal vor- 
tex core in the underdoped/overdoped materials which 
is found in this work. We also suggest possible neutron 
and ^SR experiments to probe the AF components of the 
vortex core. 

The SO (5) theory has been constructed in its general 
form to allow for anisotropy in the AF and SC couplings 
However, in the underdoped regime close to the 
AF-SC transition, most forms of anisotropies are irrele- 
vant 10]. In this work, we first study the isotropic limit of 
coupling constants, and allow only a quadratic symmetry 
breaking term. In this limit of the SO (5) theory, the free 
energy density takes the form B 
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where V' = + = 't^2 x -I- 713 y -I- 714 z are 

the SC {ip) and AF (m) order parameters. Due to the 
constraint t/'*^' + 'm^ = 1, the two are coupled. When 
g = g — Afi^X = 4x(Mc ~ M^) is negative, the bulk phase 
is superconducting; ^ > prefers the antiferromagnet. 
Assum ing a con stant direction for the Neel field, we have 
TO — a/1 — j-i/jp TO, and the equations for ip are 
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where ^ = \/p/i~g) is the coherence length; g < 
in the superconducting phase. As in the more famiUar 
Ginzburg-Landau theory of SC vortices, there are two 
length scales in the problem. 

Vortex Solutions - In searching for vortex solutions, it 
is convenient to work in polar coordinates (r, 0) and to 
rescale distance, r = ^s,and vector potential, A{r) = 



j2TT^^)s-^a{2)(i), where 



he/e* is the Lon- 



don flux quantum. The magnetic field is then B{r) = 
{(t)i^ /2t: (^)b{s) with 6(s) = s''^ da/ds. We demand 
a(0) = and a(oo) = to, the number of flux quanta 
through the plane. With ip — f (r) exp{im(j3) , then. 
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where k = Al/C is the Ginzburg-Landau parameter. We 
solve these equations by the shooting method, using the 
asymptotic solutions f{s ^ 0) ^ Cis'", a(s ^ 0) ~ 
C2S^, /(s 00) ^ 1 — C3 exp(— 2s), a{s 00) ~ 
TO — C4 exp(— 2y^/K), where Ci, . . . , C4 are constants. 
The SC order parameter and magnetic field distribution 
near the vortex core is shown in Fig. |l| for different val- 
ues of K. In the SO (5) theory, the AF order parameter 
profile is simply given by \/l — As in conventional 
vortex solutions, the SC order parameter is well approx- 
imated by a tanh. The vortex line energy has also been 
calculated and gives a value of Kcrit, separating type I 
and type II behavior, in good agreement with the value 
predicted from Hc2 for this model, Kcrit = 1 @- 

Domain Walls - The Gibbs free energy density is = 
JT— j^B-H. In the bulk, the SC state is characterized by 
I ■01 = 1, B = 0, and a free energy density of ^sc = ^ff ~ 
-p/2£_^. The AF state has ^0, B = H, and Gaf = 
-H^/Stt. Setti ug ^sc — Gaf gives the thermodynamic 
critical field H^. He — \/^TTp/£, = (I)i,/2tt^Xi^. We now 
consider a domain wall separating a bulk AF {x — > — cxi) 
from a bulk SC {x +00) and compute the energy of the 
domain wall relative to that of either bulk state. We write 
X — ^ s and A — (he/e*^) a{s) y to obtain the Ginzburg- 
Landau equations 
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FIG. 1. Vortex profiles for a series of different values of 
K. The dashed lines show the magnetic field distribution and 
the solid lines show the SC order parameter. The AF order 
parameter is related to the SC order parameter by the SO (5) 
constraint. 



with asymptotic solutions /(s —00) ^ 

Ci exp(— s^/2k), a{s —00) ^ s/k + C2, /(s — > 00) 
^ 1 — C3 exp(— 2s), and a(s 00) ~ C4 exp(— s/k). We 
again solve by the shooting method and the results are 
shown in Fig. |^ for various values of k. 

The domain wall free energy per unit length is 
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Note that the integrand of equation (||) is a difference of 
two positive quantities. When cr > the domain wall 
energy is positive. This is type I behavior. When a < 
the domain wall energy is negative and we have type II 
behavior, a — corresponds to Kcrit = 0.30, which differs 
from that determined by Hc2 or the vortex line energy 
because of the gradient which appear in the fourth order 
terms due to the S0(5) constraint. 

In the above calculations, the S0(5) constraint forces 
the vortex core to be antiferromagnetic. A normal core 
is describable within a soft superspin model. To explore 
the competition between AF and normal cores, we write 
-0 = ncos^e"^ and to = rising to with (j> = tan~^{y/x). 
We further assume that to is constant and we work in the 
extreme type-II limit where the magnetic field is ignored. 
The free energy density is then 
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where a{T) = a'{T - T^) in the vicinity of T^. We 
now consider two trial vortex profiles: (i) n(r) — 
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FIG. 2. Domain wall profiles for a series of different values 
of K. The dashed lines show the magnetic field distribution 
and the solid lines show the SC order parameter. 

no tanh(?'/^), 9{r) = (normal core), and (ii) n{r) — uq, 
6{r) — ^exp(— r/i?) (antiferromagnetic core). Here, ^ 
and £ are variational parameters, while uq = y/ (—a) /b 
is the superspin magnitude far from the vortex core. We 
find ^2 ^ 3.385 p/(-a) and = 0.9865 ^/(-g). The 
difference in free energy per unit length is then Fj^p — 
i^normai = 7rpn^(X(A) - X^), where A = 1.852y/f/a, 
= 0.3214, and 

X(A)=/ — [tanh^u-cos^de^^")] . (7) 
Jo " 

We find that the AF core is preferred for A < 0.5683, i.e. 
g/a < 0.0941. 

The consequence of this analysis is that there is a line 
in the (T, /i) plane separating regions with antiferromag- 
netic and normal cores. We find AF cores stable for 
4X(Ai2 - fi^) < 0.0941 (Tc - T), hence, AF cores should be 
observable in underdoped materials at low temperatures. 
Increasing doping or temperature will eventually result 
in normal cores. 

The general form of the free energy density when SO (5) 
has been broken down to 0(2)scXO(3)af is 

T = ip^lVV'P + ip,jVm|2 + la \rh\^ + i/3 |VP 

+ lu\m\^ + w\rn\^\^\^ + lv\^PW (8) 

where we again work in the extreme type-II limit. We 
take (3 < and u,v,w > 0. Bulk SC is stable if a/P < 
w/v; this also precludes a mixed ^ 0) phase. If 

a < and a/P > u/w we must impose > a^/u for 
global SC stability. We write ^(r) = y/\(3\/2v f {s) e"*' 
with s = {\l3\/ PjijY^^ r; f is given by the solution to 
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subject to /(O) = and /(oo) = 1. The linearized equa- 
tion for m(s) is then 

/ d\n 1 dm\ w . a 

This defines an eigenvalue problem, perhaps conveniently 
considered as a radially symmetric Schrodinger equation 
for a particle of mass M = fi p^/2p^ in an attractive 
potential V{s) — —{w/v){l — P{s))] the energy eigen- 
value is i? = {a/P — w/v). Bound states, for which 
to(oo) = 0, satisfy i? < 0, in agreement with the afore- 
mentioned conditions. Antiferromagnetic cores will exist 
for E > — T, where — T is the lowest bound state en- 
ergy; clearly T{p^/ p^,w/v) is an increasing function of 
w/v which vanishes when w/v = 0. Thus, we arrive 
at the condition w/v — T{p,^/ p^,w/v) < a/P < w/v. 
To compare with our earlier variational calculation, set 
= 1 - g/a and p^/p^ = w/v = l. 
The SC vortex with an AF core has important conse- 
quences for the high magnetic field physics in underdoped 
high Tc superconductors. Within the SO (5) theory, both 
the thermodynamic critical field = (t>-L/2n£^Xi^ and the 
upper critical field Hc2 — 4>-l/2it(^'^ describe phase transi- 
tions between SC and AF phases at fixed chemical poten- 
tial p. A schematic zero temperature phase diagram in 
the (i/, p) plane is shown is Fig. ^. Applying a uniform 
magnetic field to the AF causes the Neel vector to flop 
into the plane perpendicular to the applied field, while 
the total magnetization vector is aligned in the field direc- 
tion. The bulk AF-to-normal transition occurs at a crit- 
ical field — aJ/TipB (about 50 Tesla if a = 1), where 
J is the AF exchange constant and a is a dimensionless 
constant. Since doping (increasing p) significantly weak- 
ens J, we expect iJ^ to decrease with increasing p. On 
the overdoped SC side, we also expect i?c2 to decrease 
with increasing doping because of the loss in pairing en- 
ergy. Surprisingly, on the underdoped side, the S0(5) 
theory gives i?c2 = ^^4'l{p'^ — Pc)/'^'^P which increases 
with p. These three critical lines meet at a common tri- 
critical point Ht. Several important features are to be 
noted about our proposed phase diagram. First of all, if 
we assume that the London penetration depth Al remains 
finite at pc, then He will exceed Hc2 sufficiently close to 
Pc since behaves as p^—p^. Therefore, the SC-to-AF 
transition will change from second order to first order in 
the vicinity of ^c, where these two phases are separated 
by the thermodynamic field He cx ^fi — ^c- Secondly, all 
our discussions are carried out for a short ranged model 
at fixed p; Coulomb interactions may lead to a signif- 
icant modification of the phase diagram. However, we 
believe the most salient feature of our phase diagram, 
namely the transition from a SC state to an insulating 
state with applied field, will still remain valid in the un- 
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sign and hence generate double or multiple steps in the 
absence of other broadening mechanisms. In a longitudi- 
nal field experiment, the fluctuating transverse fields due 
to the AF vortex cores could give rise to a substantial 
amount of relaxation. The magnitude of the relaxation 
depends on the amount of spectral density of transverse 
spin fluctuations at the Larmor frequency of the muon. 
The temperature dependence of this relaxation is a probe 
of the spin fluctuations within the vortex cores. 
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FIG. 3. Schematic T = phase diagram 



derdoped regime. This could offer a basic explanation of 
the observation by Boebinger et al. ||] 

Next we consider possible experimental methods for 
observing AF vortex cores. The AF vortex core size is 
on the order of the SC coherence length (four to five lat- 
tice constants). The density of vortices is proportional 
to the field. The local electron magnetic dipolar fields in 
the cores are hundreds of Gauss, as they are in the pure 
AF |^0|, but because these AF regions are at best one- 
dimensional (along the field direction), the local fields 
may not be static, except perhaps at low temperatures. 
The correlations along these small AF tubes depend on 
the weak interlayer exchange. Spin correlations between 
AF cores of neighboring vortices will be even weaker, 
and hence the AF order may not be coherent from one 
vortex to the next. This will show up in the neutron 
diffraction pattern. For neutron scattering, one expects 
to see peaks around (tt, tt) with a width of (core size)~^. 
If the spins are not static, this scattering will be quasi- 
elastic. The total intensity can be obtained by integrat- 
ing over low frequencies. If the AF order were coherent 
from one vortex to the next, the periodicity of the vortex 
lattice would lead to superlattice diffraction peaks split 
by (inter-vortex spacing)"^. In any case, the integrated 
intensity of the broad peak at (tt, tt) will scale with field. 
For transverse field /xSR, if the AF fields are static on 
the time scale of the muon precession frequency, it may 
be possible to observe the staggered local electron dipole 
field directly. The reason is that the field distribution 
from the center of a normal vortex core appears as a step 
at the high field end of the /iSR spectrum. For AF cores, 
the dipolar field at the muon site will have a longitudi- 
nal component (along the c-axis), of order 100 G [ p^ , 
even if the electron spins lie in the a-h plane. This longi- 
tudinal field at the muon site or sites will have random 
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